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Abstract. We report results from a comprehensive study of the nearby M dwarf Proxima Centauri with the XMM-Newton 
satellite, using simultaneously its X-ray detectors and the Optical Monitor with its U band filter. We find strongly variable 
coronal X-ray emission, with flares ranging over a factor of 100 in peak flux. The low-level emission is found to be continuously 
variable on at least three time scales (a slow decay of several hours, modulation on a time scale of 1 hr, and weak flares with time 
scales of a few minutes). Several weak flares are characteristically preceded by an optical burst, compatible with predictions 
from standard solar flare models. We propose that the U band bursts are proxies for the elusive stellar non-thermal hard X-ray 
bursts suggested from solar observations. In the course of the observation, a very large X-ray flare started and was observed 
essentially in its entirety. Its peak luminosity reached 3.9 x 10 28 erg s [0.15-10 keV], and the total X-ray energy released in 
the same band is derived to be 1.5 x 10 32 ergs. This flare has for the first time allowed to measure significant density variations 
across several phases of the flare from X-ray spectroscopy of the O vn He-like triplet; we find peak densities reaching up to 
4 x 10" cnT 3 for plasma of about 1-5 MK. Abundance ratios show little variability in time, with a tendency of elements 
with a high first ionization potential to be overabundant relative to solar photospheric values. Using Fe xvn lines with different 
oscillator strengths, we do not find significant effects due to opacity during the flare, indicating that large opacity increases 
are not the rule even in extreme flares. We model the large flare in terms of an analytic 2-Ribbon flare model and find that 
the flaring loop system should have large characteristic sizes IR,) within the framework of this simplistic model. These 
results are supported by full hydrodynamic simulations. Comparing the large flare to flares of similar size occurring much more 
frequently on more active stars, we propose that the X-ray properties of active stars are a consequence of superimposed flares 
such as the example analyzed in this paper. Since larger flares produce hotter plasma, such a model also explains why, during 
episodes of low-level emission, more active stars show hotter plasma than less active stars. 
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1. Introduction 

Coronal energy release on magnetically active cool stars is, 
according to current understanding, largely driven by the in- 
terplay of surface convective motion and magnetic fields an- 
chored in the star's outer convection zone. Strong variability in 
coronal emissions due to flares (in X-rays, the extreme ultravi- 
olet, or the radio range) support theories proposing magnetic 
instabilities as the main drivers of energy release and trans- 



Send offprint requests to: M. Giidel; e-mail: 
guedel@astro . phys . ethz . ch 

* Based on observations obtained with XMM-Newton, an ESA sci- 
ence mission with instruments and contributions directly funded by 
ESA Member States and the USA (NASA) 



port (see lHaisch et aljfl99ll for a general multi-wavelength re- 
view of solar and stellar flares). Although stellar coronal radia- 
tion has traditionally been classified as "quiescent" or "flare" 
depending on the recorded time scale of variability, recent 
studies have cast some doubt on the overall meaning of this 
distinction for magnetically very active main-sequence stars. 
First, sensitivity limitations inevitably smooth out low-level 
variability. Second, a superposition of many strongly over- 
lapping flares tends to produce smooth light curves in par- 
ticular for statistical distributions of flare energies for which 
the low-energy population dominat es the total energy budget 
dKopp & Polettolll993l lGudelll2000h . And third, long time se- 
ries monitoring stellar X-ray or EUV activity have provided 
evidence for continuous variability ascribed to flares, to an 
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extent that possibly all of the observed coronal emissi on can 
be attributed to explosive energy release lAudard et alJ l2000, 
iKashvap et ai1l2002l iGUdel et alJl2003al) . This latter hypothe- 
sis is presently also in the center of solar coronal research (e.g., 
Porter et alll995tlKrucker & BenJl99llParnell & Juppfc OOO. 
Aschwanden et al 2000) n the context of "micro-" or "nano- 
"flares. 

Large, outstanding coronal flares may thus be the proto- 
types of the coronal heating agents accessible relatively easily 
by monitoring programs available on X-ray satellites. Their in- 
vestigation may provide the fundamental parameters required 
to understand explosive energy release by magnetic reconnec- 
tion, shed light on the magnetic geometry prevailing in stel- 
lar flares, and help understand processes transporting chromo- 
spheric and photospheric material into the corona. Comparing 
such studies on magnetically active stars with solar findings 
is particularly important as magnetic processes on the former 
may fundamentally differ from those on the Sun, given the 
much larger filling factor of surface magnetic fields on ac- 
tive stars. XMM-Newton and Chandra offer a new dimension 
to stellar flare research by providing access to highly sensitive 
high-resolution X-ray spectroscopy, thus potentially revealing 
information not only on the temperature structure but also on 
electron densities, turbulence, or bulk mass motions. 

During a recent XMM-Newton campaign on the nearest star 
to the Sun, Proxima Centauri, we were fortunate in recording 
a major long-duration X-ray and optical flare that outshone the 
average pre-flare emission by almost two orders of magnitude 
in X-rays, and several magnitudes in the optical. Early find- 
ings from this campaign are described in Gudel et al. (2002a = 
Paper I). The observing program was initially devised to inves- 
tigate the statistics of low-level activity, providing access to the 
weakest flares yet on any star other than the Sun. The campaign 
indeed recorded a sequence of low-level flares to an extent that 
no part of the light curve can be described as steady, for a du- 
ration of about 35 ks (Paper I). Many X-ray events were char- 
acteristically preceded by optical bursts in a manner predicted 
by the chromosphe r ic evaporation scena rio JAntonucci et alJ 
1 1984 lDennislll988l iHudson et al1ll992l) . The large flare for 
the first time provided a sufficient statistical signal-to-noise 
ratio to explicitly and significantly measure the run of elec- 
tron densities in the flare plasma from spectroscopic He-like 
line triplets. Electron densities of several times 10 11 cirT 3 were 
found during peak time, while they rapidly decreased during 
the flare decay phase. These observations compare favorably 
with previous solar spectroscopic measurements of flare densi- 
ties jMcKenzie et alJl980llDoschek et alJl98l . 

The present paper presents the complete XMM-Newton data 
set on Proxima Centauri and disc usses some models for the 
large flare. In a companion paper llReale et alJl2004l) . we will 
present a series of hydrodynamic simulations for the same flare. 
The paper is structured as follows: Sections 2 and 3 describe 
the target and the observations, Sect. 4 presents the analysis 
and our results, while Sect. 5 discusses the findings in the con- 
text of flare physics and coronal heating. Sect. 6 contains our 
conclusions. 



2. The Target 

Proxima Centauri is a magnetically active dM5.5e dwarf (V = 
11.05, see SIMBA D database; L bol = 6.7 x 10 30 ergs" 1 , see 
iFro gel etalJfl972h revealing strong coronal activity. At a dis- 
tance of 1 .295+0.005 pc JPerrvman et all 19971) . it is the closest 
star to the Sun. Its "quiescent" X-ray luminosity Lx ~ (4- 16)x 
10 26 erg s' 1 faaisch etalJl990l) is similar to the Sun's despite 
its «50 times smaller surface area. It h as attracted the atten tion 
of most prev i ous X -ray observatories. Hai sch et al.l IU ,980|) and 
lHaisch et alJdl98ll) observed quiescent and flaring X-ray emis- 
sion from this star with the Einstein satellite, with the largest 
flare reaching an X-ray luminosity of Lx = 7.4 x 10 27 erg s _I 
and peak temperatures of 17 MK during the rise to the luminos- 
ity peak. The surface filling factor of X-ray emitting magneti- 
cally confined plasma was estimated to be about 20% based on 
solar analogy. Interestingly, coordinated radio, optical, and ul- 
traviolet observations recorded no flare emissions. The authors 
interpreted this long-decay flare (decay time scale t * 20 min) 
by a series of arches cooling predominantly by radiation, with 
dimensions compa rable to the stellar radius. In a follow-up ob- 
servation, Hai sch et alJ (119831) report another very large flare, 
with a peak X-ray luminosity of Lx = 1.4 x 10 28 erg s _1 in 
the 0.2^4 keV band, and a maximum temperature of 27 MK. A 
two-ribbon flare model in analogy to large solar flares was pro- 
po sed for this long-du ration event (r =s 20 min), refined later 
bv lPoletto et alJ {l]?88). Further modeling of this flare, involv- 
ing hydrodynamic simulations of flaring loops, were presented 
bv lReale e t al. ( 1988). Both latter references found flaring loop 
structure s with s i zes up to about 7 x 10 9 cm or 0.7 R*. 

lHaisch et al] Jl990l) present a coordinated EXOSAT-IUE 
observation of Proxima Centauri, again comprising o ne larger 
flare (peak L x = 3 x 10 27 erg s 1 ). Hai sch et alJdl995l) obtained 
rather sensitive observations of Proxima Centauri in the harder 
bands provided by the ASCA detectors where the contrast be- 
tween (hotter) flares and the (cooler) quasi-steady emission is 
more pronounced. The authors point out that a large number of 
intermediate flares seem to be present in the light curve, corre- 
sponding to solar M-class flares. Even harder emission was de- 
tected by the XTE satellite, but unfortunately no large flares oc- 
curred during the time of the observations (Hais ch et alJl9 98). 

Proxima Centauri is an obvious target also for the newest 
generation of satellites, XMM-Newton and Chandra. First re- 
sults of the observation present ed here have been reported in 
Paper I and lGiidel et all J2003b) pertaining to flare density in- 
creases, the Neupert effect (Neupert 1968) for X-ray and opti- 
cal fla res, and continuous low-level flaring. Iwargelin & Drake! 
(2002) discuss a CCD observation obtained by Chandra con- 
centrating on the possibility to detect charge exchange from the 
stellar wind. 

3. Observations 

XMM-Newton Jjansen et alJl200ll) observed Proxima Centauri 
on 2001 August 12 during 65 ks of exposure time. All detec- 
tors were operating nominally. The observing log is given in 
Tabled The data were reduced using the Stan dard Analysis 
System (SAS), version 5.3 for the EPIC data dStruder et alJ 
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l200lllTurner et all200ll) and version 5.4. 1 for the RGS spectra 
Jden Herder et all200ll) . 

In anticipation of possible large flares occurring on this star, 
we used the small window modes for all EPIC detectors, avoid- 
ing X-ray pile-up effects up to count rates of about 5 cts s _1 for 
the two MOS detectors and up to 130 cts s -1 for the PN camera. 
The large flare reported here exceeds these limits for the MOS 
but not for the PN. We will thus use only the PN for CCD spec- 
tral analysis but combine data from all three detectors to pro- 
duce sensitive light curves during the low-level emission. 

We extracted the EPIC spectral data from circles with a ra- 
dius of 50" (encircling 90% of the source photons). For the 
PN, the corresponding response matrix and the effective area 
were derived using the standard SAS tasks rmfgen and arfgen. 
The energy range of sensitivity of the EPIC cameras is 0.1- 
12 keV. Their energy resolving power is best at highest ener- 
gies; it reaches E/AE * 50 at 7 keV for the PN and scales as 
E 1 ' 2 . 

Background light curves were extracted from CCD fields 
outside the source region (for the MOS, we used four back- 
ground areas in the outer chips for improved statistics; back- 
ground vignetting is negligible since the background ist domi- 
nated by intrinsic detector noise and by soft proton irradiation; 
both components are only weakly affected by vignetting; we 
refer to the XMM-Newton Users Handbook, § 3.3.7 for fur- 
ther details). A few small background flares induced by pro- 
ton events in the CCD were seen, with amplitudes in the com- 
bined EPIC PN+MOS1+MOS2 data (see below) of no more 
than 0.65 ct s _1 during the first 1.7 hrs, and of no more than 
0.33 ct s _1 thereafter. These events were reliably eliminated 
by subtracting the background light curve. The median back- 
ground rate was relatively high, with 0.08 cts s _1 in the com- 
bined light curve. 

The two Reflection Grating Spectrometers (RGS) provide 
high-resolution spectroscopy (E/AE « 300 FWHM around the 
Ovn lines at « 22 A) over the wavelength range of 5-35 A, 
with a maximum first order effective area of about 120 cm 2 at 
15 A. To include recent updates of the effective area function, 
we reanalyzed the RGS data using the SAS 5.4.1 release but 
note that there is no significant change with respect to the re- 
sults reported previously (Paper I, using SAS vers. 5.3.3 for the 
data reduction). The background was determined from adja- 
cent areas on the RGS CCD detectors as prescribed in the SAS 
task rgsproc. The response matrix (including the effective area 
function) was produced with the SAS task rgsrmfgen. 



Table 1. Observing log 



2001 August 12 



Instrument 


UT ranj 


;e 








JD 2452133.5+ 


PN 


04 : 


: 21 


: 29 


-22 : 


: 37 : 


: 29 


0.18147 


- 0.94270 


MOS1 


04 


: 01 


: 29 


-22 : 


: 35 : 


: 08 


0.16770 


- 0.94106 


MOS2 


04 : 


: 01 


: 29 


-22 : 


: 35 : 


: 08 


0.16770 


- 0.94106 


RGS1 


03 


: 53 


: 59 


-22 : 


: 37 : 


: 34 


0.16249 


- 0.94275 


RGS2 


03 : 


: 53 


: 59 


-22 : 


: 37 : 


: 32 


0.16249 


- 0.94273 


OM 


03 : 


: 58 


: 58 


- 17 : 


: 43 : 


: 46 


0.16595 


- 0.73873 



The Optical Monitor (OM; Mason et al.l200lh observed the 
star through its U band filter in high-time resolution mode. The 
data reduction used the updates provided by SAS vers. 5.3.3. 
The count rates observed during the flares approach saturation 
of the detector. As the true photon flux must be estimated based 
on detector dead-time corrections, the photon count rate num- 
bers given in this paper are somewhat uncertain and are based 
on the current understanding of the calibration. Given the un- 
certain physical nature of the optical emission, the level of ac- 
curacy is entirely adequate for our purposes. 

4. Data analysis and results 

4.1. Light curves 

The light curves extracted from the EPIC cameras are shown in 
Fig. 1 and 2. The first «T2 hours of the observing time recorded 
low level radiation. This portion of the observation is shown in 
Fig. 1. At t sa 0.715 d, an extremely large flare sets in, rising to 
its luminosity peak in about 12 minutes. A very long decay fol- 
lows, with several superimposed secondary flares visible. This 
episode is illustrated in Fig. 2. 

For our later analysis, we define time sections of i) pre-flare 
low level emission (Q) ii) flare rise (A), iii) first flare peak (B), 
iii) initial flare decay (C), iv) secondary flare peak (D), and iv) 
final flare decay (E), as shown in Figs.n]and|2] 

Fig. 1 shows, from top to bottom, the summed EPIC 
PN+MOS1+MOS2 count rate at a resolution of 60 s. We have 
confined our photon extraction to the energy intervals 0.15- 
4.5 keV for PN and 0.18-4.5 keV for MOS in order to sup- 
press contaminating warm pixel contributions and excessive 
background, while optimizing the signal-to-noise ratio. Tests 
with lower energy limits up to 0.4 keV did not show signifi- 
cantly different results. The initial ^20 min (until 0.1815 d, see 
Table Q were not observed by the PN camera; we multiplied 
the MOS count rates by 2.7 to obtain an estimate for the total 
PN+MOS light curve for that interval but note that, due to the 
spectral dependence of the PN/MOS count rate ratio, this esti- 
mate is only approximative and the absolute count rate values 
before 0.1815 d should be accepted with some reservations. We 
therefore also plot the MOS-only light curve for the early time 
interval (dashed in Fig.Q. The middle panel in Fig. 1 presents 
the hardness ratio at a resolution of 300 s, here defined as the ra- 
tio between the count rates in the 1 — 4.5 keV and the 0.4-1 keV 
bands. Hardness is thus a sensitive indicator of heating events. 
The lowest panel shows the OM light curve at a resolution of 
80 s per bin, chosen for optimum illustration of rapid flares. 

The most surprising aspect of this part of the observation 
is the level of continuous variability. At no time is the light 
curve constant; even if the long-term trend is removed, there 
are hardly any intervals of more than a few tens of minutes 
that are free of modulations at the given sensitivity. This is pro- 
viding evidence of continuous microvariability. There are three 
obvious time scales of variability: 

i) The X-ray light curve slowly decays by an order of 
magnitude from an initial ss 10 cts s to w 1.5 cts s at 
t = 0.7 d. The count rate-to-luminosity conversion factor ap- 
plicable to this part of the light curve has been derived from 
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0.20 0.30 0.40 0.50 0.60 0.70 

UTtime on 2001/08/12 (d) 



Fig. 1. Low level emission of Proxima Centauri. Top panel: X-ray light curve. Counts from all three EPIC detectors have been 
co-added (using counts between 0.18^.5 keV from the MOS detectors, and 0. 15-4.5 keV from the PN). Time resolution is 60 s. 
The horizontal bar marks the interval that was used for spectral analysis of low-level emission. The vertical lines (also in the other 
panels) indicate flare-like features. Solid lines refer to clearly detected features (in the X-ray or the optical band), while dotted 
lines mark borderline cases. The dashed curve until 0.22 d is the MOS-only light curve shown here because the PN was not in 
operation before « 0.18 d (the total light curve has been approximately scaled). Middle panel: Hardness, defined as the ratio of 
counts between 1-4.5 keV to counts between 0.4-1 keV, at a time resolution of 300 s. Again, counts from all EPIC detectors 
were co-added. Bottom panel: Optical U-band light curve from OM (on a logarithmic count rate scale). The time resolution is 
80 sec per bin. Time gaps are of instrumental origin, specific for the selected observing mode. 



spectral modeling and amounts to w 2.8 x 10 erg s per 
ct s _1 . From this, the low-level luminosity varies between « 
(4 - 28) x 10 26 erg s _1 , compatible with previous literature val- 
ues. The cause of the long-term decay is unclear. It could be 



due to rotational modulation of a stellar surface that is inhomo- 
geneously covered by X-ray emitting active regions, most of 
which rotate out of sight during the observation. The rotation 
period of Proxima Centauri is, however, suspected to be as long 
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Fig. 2. Similar to Fig.^ with the following exceptions: The top and the middle panels additionally show the PN-only curves 
(dashed), since only those data were not piled up during the flare peak. And the optical light curve in the bottom panel is shown 
at a time resolution of 40 s, to illustrate some finer details during high flux. The OM observation was terminated early, but covers 
the main flare peak. The intervals A-E used for spectroscopy in this paper are marked in the top panel. 



as 83.5 d, see lBenedict et alJ dl998l) . Alternatively, the slowly 
decaying emission may be part of the late light curve of a giant 
flare occurring mostly before our observations, similar to the 
flare occurring during the second part of our observations. 

ii) A shorter time scale of about 1 hr reveals itself as a 
slow modulation in the X-ray light curve. At closer inspec- 
tion, some of gradual peaks are in fact composed of superim- 
posed weak flares, as was shown in Paper I. They often cor- 
respond to groups of simultaneously occurring optical flares. 



A particularly clear case is seen at t — 0.43 - 0.46 d in Fig. 
1. The hardness light curve mostly shows enhancements dur- 
ing the increasing portions of the X-ray light curve on the same 
time scale. This indicates increased flare activity during the rise 
phase and decreased activity during the decaying phase when 
plasma cools and the emission becomes softer. 

iii) We eliminated the above modulations from the light 
curve by subtracting a smoothed version of the curve, leaving 
only short-term modulations in the residuals. The autocorrela- 



tion function of this "high-pass" filtered curve shows a decay 
from unity at zero lag down to zero correlation at about 4 min- 
utes time lag. This is not compatible with residual noise but 
in fact corresponds to the typical time scales of the weakest 
flares, some of which were illustrated in Paper I. These flares 
have peak luminosities of (2 - 5) x 10 26 erg s _1 , corresponding 
to total released soft X-ray energies of « 10 28 - 10 29 erg (in the 
0.1-10 keV range). Such flares would be classified as "M-class 
flares" on the Sun (see lHaisch et al]l 19951 for a discussion of 
M-class flares in the context of stellar observations). Many of 
the detected X-ray flares are characteristically accompanied by 
optical flares which typically occur during the rise time of the 
X-ray flare and provide support for a flare model in which ac- 
celerated electrons (signified by the optica l burst as they collide 
in the chromosphere; Huds on et alJI 19921) deposit their energy 
in the chromosphere, heating and evaporating t he plasma that 
subsequently fills the magnetic loops (see, e.g. jAntonucci et all 
11984 lDennisi ri988). Such "chromospheric evaporation" ap- 
pears to occur at a high cadence in Proxima Centauri; clear 
examples are marked by solid vertical bars in all three panels 
in Fig. 1, somewhat less obvious cases by dotted bars. Most of 
these events have been described in Paper I; we refer to that pa- 
per for a more detailed discussion on low-level flaring. Some 
of the outstanding examples can also be traced in the hardness 
curve although the dominating continuous or slowly varying 
level of radiation overwhelms short-term hardness changes in 
the superimposed flares. 

Fig. 2 shows the large flare, in three panels equivalent to 
those described in Fig. 1, except that we plot the optical light 
curve with a time resolution of 40 s to bring out some finer de- 
tails during the high-count rate episode. The Proxima Centauri 
source was partly piled up in the MOS detectors during the 
flare, reaching »23 cts s _I each at flare peak, while the upper 
limit to avoid CCD pileup is « 5 cts s _I for the MOS small win- 
dow mode (see XMM-Newton Users' Handbook, Sect. 3.3). 
The pile-up limit for the PN small window is 130 ct s _1 , while 
our flare peaks at * 78 ct s _1 . We therefore also plot the 
count rate and hardness curves for data extracted from the PN 
detector only. The total count rate and the hardness curves 
appear to be little affected by MOS pileup, however, as the 
comparison between the light curve pairs shows. The appro- 
priate count rate-to-luminosity conversion factor for the PN- 
only data is 4.95 x 10 26 erg s per ct s _1 at flare peak. The 
ratio changes only within a few percent from the flare peak 
to the decay phase. The peak luminosity of this flare is thus 
ss 3.9 x 10 28 erg s _1 [0.15-10 keV], and the total energy emit- 
ted in X-rays is of order 1.5 x 10 32 erg. As predicted by the 
chromospheric evaporation scenario, hardness peaks during the 
rise of the X-ray light curve, i.e., the highest temperatures are 
attained before the emission measure has built up to its maxi- 
mum. Similarly, the optical light curve peaks during the most 
rapid increase of the X-ray curve. We also note that both the 
X-ray and the optical light curves show a very faint precursor 
flare, marked by a vertical line, and again it is the optical signal 
that precedes the X-ray flare. The slowly decaying light curve 
shows two large, superimposed events, one at 0.75 d and one at 
0.85 d. Unfortunately, the OM instrument did not continue to 
observe sufficiently long to cover these two events as its moni- 
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Fig. 3. EPIC PN spectra during various stages of the flare (A-E) 
and during the pre-flare episode Q. The count rates have been 
multiplied by various factors to prevent mutual overlap (see 
factors indicated at lower left). Note the variable high-energy 
slope, indicative of variable flare temperatures. 



toring program was terminated prematurely for operational rea- 
sons. They differ in that the earlier one does not result in a 
significant increase of the hardness (but rather in a plateau of 
constant hardness), while the later flare does. This may be the 
result of "hardness contrast": When the later secondary flare 
occurs, the hardness of the decaying initial flare plasma has 
already decayed nearly to levels seen in Fig. 1. Further faint 
flares can be seen during the decay phase, such as the one at 
0.91 d. They seem to correspond to ongoing low-level flaring 
as detected in Fig. 1. 

4.2. Spectroscopy with EPIC and RGS 

The EPIC PN spectra shown in Fig. [3] characterize the spec- 
tral evolution of the large flare across episodes A-E (together 
with the low-level spectrum 'Q'). The individual spectra have 
been shifted along the y axis by factors indicated in the figure 
to avoid mutual overlap. The high-energy slope (between 1.5- 
7 keV) is largely determined by the bremsstrahlung continuum 
of the hottest dominant plasma component. As the flare pro- 
gresses, the slope becomes steeper, indicative of overall cool- 
ing of the flare plasma. Also, a prominent Fe K complex in- 
cluding He-like Fe xxv lines is seen around 6.7 keV predomi- 
nantly in the earlier flare phases and at flare peak. Again, these 
lines are indicative of the presence of very hot (> 10 MK) 
plasma. And third, the Fe L-shell region around 0.6-1.0 keV 
is another characteristic indicator of temperature since it con- 
tains Fe lines with ionization stages between Fe xvn-Fe xxiv 
that shift the peak emission of the unresolved spectral maxi- 
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mum from « 0.75 keV in the pre-flare spectrum to «1 keV at 
flare peak, decreasing again as the flare progresses. Various line 
features above 1 .4 keV can be used to constrain abundances of 
Mg, Si, S, Ar and Ca. We will present an abundance analysis 
below, combining EPIC information with RGS spectroscopy. 

Reflection grating spectra obtained during the low-level 
episode before the flare and during four intervals from flare 
peak to late decay are shown in Fig.|4] (The spectrum obtained 
from interval A shows a very low signal-to-noise ratio; the anal- 
ysis pertaining to this interval strongly relies on the EPIC PN 
spectrum.) A few important lines are marked. The spectra ob- 
tained during low-level emission (interval Q) and at flare peak 
(interval B) are strikingly different, indicating heating to high 
temperatures during the flare. For example, a marked enhance- 
ment of the continuum flux level is visible at flare peak, and 
the appearance of lines of Fe xx-Fe xxiv around 10-12 A with 
peak formation temperatures between 10-20 MK indicate the 
presence of large amounts of hot plasma. 

The spectra were used to study emission measure (EM) dis- 
tributions, abundances, densities, and optical depths at various 
stages of the flare development. We discuss results in the fol- 
lowing subsections. 

4.3. Spectroscopic density measurements 

The He-like line triplets of Ovn (at 21.6-22.1 A) and of 
Neix (at 13.4-13.7 A) can be used to obtain characteristic 
electron densities in the source region of the respective lines 
(Gab riel & Jordan! 19691) . The ratio of the fluxes of the forbid- 
den line (/, 1 s 2 3 S o - 1 s2s 3 S i , at 22. 1 A for O vn and at 1 3.7 A 
forNeix) and the intercombination line (i, Is 2 1 Sq- ls2p 3 P\ 
at 21.8 A for O vn and at 13.55 A for Neix) are sensitive to the 
electron density n e in the range of 10 9 - 10 12 cm" 3 for Ovn 
and 10 10 - 10 13 cm" 3 forNeix IPorauet e t~alll200ll) . Further 
He-like triplets, in particular of Mg xi and Si xm, are accessible 
within the RGS range, but the spectral resolution is insufficient 
to cleanly separate the the lines of the triplet, and the density- 
sensitive ranges (> 5 x 10 1 1 crrT 3 ) are probably only marginally 
useful for typical coronal flare plasmas on dwarf stars. 

In Paper I, we discussed the evolution of the f/i ratios dur- 
ing the large flare episode; we will not repeat the analysis but 
briefly summarize the results. The densities derived from Ovn 
flux ratios were found to peak around 4 x 10 11 cirT 3 during the 
primary and secondary flare peaks (intervals B and D, respec- 
tively), while n e was lower by about one order of magnitude 
during the decay episodes (intervals C and E). These densi- 
ties combined with the spectral measurement of the EMs of the 
O vn line forming plasma were used to estimate total masses 
and volumes during the flare. It was found that the mass of the 
relatively cool O vn forming plasma increased appreciably to- 
ward the decay episodes, possibly due to continuous cooling of 
plasma that was initially heated to much higher temperatures. 
The O vn line triplets are represented in flux units in Fig.[5](B- 
E). Suppression of the / line in panels B and D is clearly vis- 
ible. The spectrum obtained during the low-level episode (Q) 
shows a line ratio compatible with fluxes around the lower limit 



of the density-sensitive range of the O vn triplet, corresponding 
to densities in the range 10 9 - 10 10 cm -3 (see Paper I). 

The analysis of the Neix triplet (see Fig.[6]for an illustra- 
tion of fluxed spectral extractions) is much more complicated 
due to stro ng blending with several lines in particular of Fe xix 
jNess et al . 2003b, Fe xix has a peak formation temperature of 
about 8 MK, compared to 4 MK for Ne ix). While Fig. |6jT and 
E do show a flux peak at the position of the / line at 13.7 A, its 
possible suppression in panels B and D cannot be established 
given the dominance of contributions from Fe xix at flare peak. 

In principle, there may be some contributions to the i and 
/ lines from the steady component; the // i ratios may thus be 
slightly biased by non-flare contributions. Since the / line is 
much stronger than the i line during the low-level emission 
(Fig. 02), the correction is essentially for the former flux. If 
the decaying trend of the overall light curve (Fig.^ continues 
after the start of the large flare, then, for the intervals D-E, we 
expect an Ovn flux about 2.5 times smaller than on average 
during the low-level interval Q, resulting in a non-flare / flux 
about 13% of that observed during interval D. Correction of 
this effect would increase the density by about 15%, which is 
not significant given the considerable errors in the flux ratios. 
In any event, we note that the densities reported in Paper I are in 
this sense lower limits, but specific corrections are meaningless 
since we cannot accurately estimate the (small) contribution 
of the non-flare emission at any time during our strong flare. 
The effect of low-level contributions of course becomes more 
significant at later times during the flare decay. We have also 
investigated the f/i behavior during the late secondary flare af- 
ter 0.84 d, but no significant density increase as in the earlier 
episodes could be observed. 

It is interesting to note that the density behavior of this 
flare finds para l lels in observations of large solar flares. 
iMcKenzie etaD ([1980) presented Ovn flux ratios during a 
large flare with similar time scales as ours and found f/i ra- 
tios around unity close to the flare peak, implying densities 
of up to 2 x 10 11 cm" 3 . Much shorter flares were discussed 
bv lDoschek et alJ dl98lh : in those cases, the densities reached 
peaks around (10 - 20) x 10 11 cirT 3 as measured from Ovn, 
but this occurred very early in these flares (during the flare 
rise), while during flare peak, electron densities of a few times 
10 11 cirT 3 were derived. As in our observations, the densi- 
ties very rapidly decreased after the peak to levels at or below 
10 11 cm" 3 . The estimated masses and volumes, on the other 
hand, steadily increased, a conclusion also drawn in Paper I for 
Proxima Centauri where we suggested that this indicates an in- 
crease of the amount of cool plasma that was initially heated to 
higher temperatures but that cools to a few MK. Unfortunately, 
we have no reliab le density diag nostic for the more relevant 
flare temperatures. lDoschek et alTdl98ll) and references therein 
report solar-flare densities derived from Fexxv (T > 10 MK) 
that are similar to those measured from Ovn (T * 2 MK). 
lEandietal] §003) recently used various density diagnostics 
for an M class solar limb flare. From Fe xxi lines, they derive 
densities up to 3 x 10 12 cm" 3 , but there are conflicting mea- 
surements for lower ionization stages that reveal much lower 
densities, comparable with pre-flare densities (see also refer- 
ences in that paper). If the density values at « 10 7 K are real, 
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wavelength (A) 



Fig. 4. Fluxed RGS1+2 spectra referring to the intervals Q (pre-fiare low level), and B-E (see Fig. [2] for reference). The bin 
resolution is 50 mA except for panel B where a resolution of 78 mA has been used to improve the signal-to-noise ratio. 
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then pressure equilibrium cannot be assumed for flaring loops; 
a possible explanation involves volumes for the O vn and the 
Fe xxi/Fe xxv emitting plasmas that are spatially separate, pos- 
sibly contained by different magnetic field lines. Obviously, 
further solar studies are needed. 

4.4. Abundances and emission measure distributions 

A derivation of elemental abundances is inherently related to 
the knowledge of the EM distribution because the formation 
efficiency of all observed X-ray emission lines in a collision- 
ally ionized plasma is strongly dependent on the tempera- 
ture structure. For previo us efforts using high-resolution X- 



(2001j.lGudel et alJl2001ak andlAudard et all 


(2003a). EM dis- 


tributions are also discussed in 


Huenemoerder et alJ 


(2001). 



Despite some differences between the reported results, there 
is wide agreement that magnetically very acti ve stars reveal 
an in verse First-Ionization Potential effect (Hrinkman et alJ 
l200ll) that expresses itself in low abundances of elements 
with a low first ionization potential (FIP), and enhanced abun- 
danc es of elements with a high FIP. This is in contrast to the 
Sun dvon Steiger & Geissl il989) and solar-like inactive stars in 
which a normal FIP effect has been found i n similar X -ray stud- 
ies (enhanced low-FIP elements; Giidel et al. 2002bJ). 

There is particular interest in studying elemental abun- 
dances also during flares. Flares bring chromospheric and pho- 
tospheric material into the corona after heating this cool gas 



and inducing chromospheric evaporation. Differences between 
the composition of flare plasma, chromospheric plasma, and 
the overall low-level coronal plasma could have important im- 
plications for mass transport or diffusion processes. While 
overall abu ndance changes in flares have been reported previ- 
ously (e.g., iMewe etalJll997h . a selective FlP-related change 
of elemental abundances was first reported for a large flare 
on an RS CVn-type binary bv lGiideletal.l(ll999l) . Similar en- 
hancements of low-FIP elements duri ng large flares on other 
RS CVn binaries w ere described by Ost en et alJ (|2000) and 
lAudard et alJ J200ll) . We will test here whether similar en- 
hancements are seen in the Proxima Centauri flare. 

Despite the comparatively large count rates in the present 
flare on Proxima Centauri, the analysis of abundance anoma- 
lies and the differential EM distribution is challenging because 
of the short time scales of variability. We again use sections 
Q and A-E marked in Fig.[0for the analysis. We subjected the 
PN and R GS 1+2 data to the same analysis method as described 
in lAudard et alJ (12003 at) . In short, we used only the strongest 
emission lines of the RGS and their immediate surrounding 
wavelength ranges, in particular lines of Mg, Fe, Ne, O, N, and 
C, plus a few intervals of continuum. We thus preserve a correct 
treatment of the overlapping line wings in particular in the Fe 
L-shell region while discarding many regions containing weak 
and ill-defined emission lines with often poor atomic data. We 
also flagged wavelength intervals with poor instrumental cal- 
ibration, in particular all data shortward of about 8.3 A in the 
RGS. The description of the important hot flare plasma compo- 
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nents requires additional information from the EPIC PN cam- 
era that sensitively constrains the hot end of the EM distribu- 
tion by the bremsstrahlung continuum and the Fe K complex at 
6.7 keV. Since the high-resolution RGS data discriminate much 
better between spectral features at lower energies, we used the 
PN data only above 1.3 keV (shortward of 9.3 A). The com- 
bined PN + RGS1 + RGS2 data sets were the n subject to a 
multi-c omponent analysis in the SPEX software dKaastra et alJ 
1996a) to derive an initial set of temperatures, EMs, and abun- 
dances. This procedure has the advantage of optimally consid- 
ering all line blends and wings that may affect the selected line 
systems, while at the same time self-consistently calculating 
the (abundance and temperature-dependent) continuum level. 
The abundances were then used to derive a more comprehen- 
sive EM distribution. In principle, this procedure could be it- 
erated, but the limited data quality at hand does not warrant a 
deeper analysis. 

We report here the abundance ratios (with respect to the Fe 
abundance) and EM distributions obtained from the intervals 
Q and A-E. All ratios are normalized to sol ar photospheric ra- 
tios as given by Ander s & Grevessd (LL989). For Fe alone, we 
adop ted an updated solar p hotospheric abundance reported by 
iGrevesse & Sauval dl999l) . We prefer to use these traditional 
solar values for ease of comparison with previous X-ray stud- 
ies. We do, however, also show the abundance ratios if the 
upd ated solar photospheric value s of C and O presented by 
lAllende Prieto et alJ d200ll Eo02) are used as a basis. The re- 
sults relative to the former set of solar photospheric abundances 
are listed in Table[2]together with their lcr rms errors, and illus- 
trated in Fig.0(filled circles). Since the updated Allende Prieto 
et al. photospheric abundances of C and O are smaller by fac- 
tors of 1 .48 and 1 .74, respectively, than those given by Anders 
& Grevesse, the ratios in Fig.0increase by the respective fac- 
tors. These modified abundance ratios are shown by open cir- 
cles. 

The most surprising aspect here is that the abundance ratios 
with respect to Fe are nearly constant across the range of FIP, 
and most ratios are actually close to solar photospheric values 
(ratios of unity). The abundances of Ne, N, and C tend to sys- 
tematically exceed unity, although the effect is, given the error 
bars and unknown systematic errors in the atomic physics (in 
particular of several blending lines), rather small . If the upd ated 
solar abunda nces for C and O given by Allende Prieto et al. 

are considered, then there is some trend during 
the late flare phases for high-FIP elements to be higher than 
low-FIP elements (by factors of 2-3), while this trend is not 
visible during flare peak. We also quote the "absolute" abun- 
dance (with regard to H) of Fe for each interval; such values 
require good measurements of the continuum which we feel 
confident we have well determined using the high-energy tail 
of the EPIC PN spectrum (see Fig.|3|l, given the high dominant 
temperatures in the flare (Figs.|8j|9j. 

Fig.[8]and[9]illustrate two versions of the corresponding EM 
distributions. The former set was obtained by fitting Chebychev 
polynomials of order 6, while the latter set uses an integral 
inversion technique with a regularization parameter for the 
smoothness of the DEM (errors cannot be derived in the con- 
text of the polynomial fit, but comparing with the alternative 
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Fig. 7. Elemental abundances determined during the low-level 
episode Q and during flare intervals A-E. All abundances 
are plotted as ratios to the Fe abundance, and are normal- 
ized with th e respective solar photos pheric abundance ratios 
according to Anders & Grevesse ( 1989). For Fe, an updated 
value given bv lGrevesse & Sauvall 41999) has been used. The 
op en circles refer to solar a bundanc es of O and C reporte d 
bv lAllende Prieto et alJ J200lh and lAllende Prieto et alJ hOO^t . 
respectively. 



regularization method illustrates the range of solutions). The 
differences between the methods are not surprising and relate 
to the ill-posedness of the underlying integral i nversion prob- 
lem. De tails of these methods are described in iKaastra et alJ 
( 1996b). A comparative study of several inversion techniques is 
given in Giide l et all {l297a). Overall, the polynomial method 
tends to reveal peaks in the EM distribution, while the regular- 
ization method shows extended wings up to high temperatures. 
The two versions illustrate the principal uncertainty in the re- 
construction of an EM distribution but at the same time clearly 
illustrate the main thermal structure of the flare, and its change 
across the different flare phases. The lowest panel in each of 
these figures illustrates the evolution of the EM distribution by 
showing all DEMs on the same scale. Although the regulariza- 
tion method clearly spreads out the EM over a larger range than 
the polynomial method, the trends are the same: Rapid heating 
to high temperatures is followed by an increase of the EM to 
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Table 2. Abundance ratios with respect to Fe for various phases (A-E) of the large flare (values rounded to multiples of 
0. 1), and for the low-level pr e-flare section (Q). All values are normalized with the sola r ratios. The basic solar a bundances 
of lAnders & Grevessd (l989) have been used, except for Fe for which the value given in iGrevesse & Sauval i 19991) has been 
adopted. Errors are based on lcr errors of the individual abundances (see text for details). "Fe/Fe" is unity by definition, but 
the error bars indicate the fractional uncertainty of the Fe abundance. The last line gives the Fe abundance relative to solar 
photospheric values. 



Element 


Q 


A 


B 


C 


D 


E 


C/Fe 

N/Fe 

O/Fe 

Ne/Fe 

Mg/Fe 

Si/Fe 

S/Fe 

Ar/Fe 

Ca/Fe 


1 n+(>7 

1 5 +a6 

1,J -0.6 

6 +0 2 

"•"-0.2 

1 6 +06 

7 i +0.9 

"^-o.s 

7 i +0.8 
A1 -0.8 


L8 -L9 

1 7 + " 

1 r - 6 

1 ' 1 -0.7 

2 4 +2 - 2 

-2 1 

2 +0 - 9 

"•^-0.2 

2 +04 

"•^-0.2 

6 +0 - 6 

"•"-0.5 


°- 7 -03 

i r 4 

5 +a2 
"■ J -o.i 

1 3 +0 - 4 

6 +0 2 

-0 2 

1 +02 

1 -"-0.2 

7 +a2 

"• ' -0.2 
"•"-0.3 


L "-04 
i 9+0.4 
I - z '-0.4 

1 +0 2 

I -"-0.2 
L8 -0.4 

1 2 +0 - 3 
1.2^ 

o.8!iJ 2 

5 +0 5 

"• J -0.4 


I -7+0 4 

II '-0.4 

1 5 +0 - 4 

i r - 2 

1 ' 1 -0.2 

1 6 +03 

l.u 3 

1 6 +0 - 3 
1 5 +0 - 2 

lp -0.3 

i r - 2 

1 - 1 -0.2 

9 "3+0.7 
^ -0.7 

1 5 +L0 


n+0 4 

2.0 + _Z 

1 q+0.4 
i-^-0.4 

1 2 +0 - 2 

1 - z '-0.2 
9 q+0.4 
^••'-0.4 

1 3 +0 - 3 

1- -0.3 

1 3 +0 - 2 

lp -0.2 

i o +0 - 3 

l .o o .3 
"• J -0.3 


Fe/Fe 
Ni/Fe 


o.r- 


1 +0 ' 4 


1 o +0 - 2 

i -"-0.2 
q i +0.9 
3 -0.9 


1 +0 ' 2 

i -"-0.2 

3 6 +u 
•'•"-l.o 


1 +01 

-0.1 

4 r 10 

^• x -i.o 


1 +01 

9 7+O.8 
^•'-0.8 


Fe/H 


0.51± 


0.58!S 


, 91+0.20 


"■"'-o.io 


54+ 007 

J -0.06 


"•^ -0.05 



its peak at somewhat lower temperatures, after which the EM 
spreads out over a larger range while cooling. 

4.5. Optical depth 

Standard interpretation of coronal emission assumes that the 
optical depth of the plasma is negligible for both line and con- 
tinuum emission. Traditional spectral analysis of coronal EM 
distributions and elemental abundances would be fundamen- 
tally flawed if significant optical depths were p resent, and ra- 
diative transport methods would be in order. ISchriiver et alJ 
(1994) argued for non-zero optical depth due to resonant scat- 
tering in several lines of the inactive solar analog a Cen in 
the extreme ultraviolet range, a view that was challenged by 
a fol low-up investigation including EM analysis from X-ray 
data ( Schmitt et alJll996l). The situation is also rat her unclear 
in the solar context. While lSchmelz et alld 19971) and lSaba et alJ 
i 19991) find evidence for optical depth effects in the Fexvn 
/115.01 line based on flux ratios with respect to the Fexvn 
/115.26 and the Fe xvn /116.78 lines, recent laboratory measure- 
ments of these r atios sign ificantly differ fr om previous fheoreti - 
cal calculations jBrown et alJl998lEoOlllLaming et aJ20 00). 

Stellar coronal optical depths have recently been stu died 
using grating observati ons from XMM-Newton jAudard et alJ 
l2003al) and Chandra JNess et al.ll200lL l2003al) . with results 
compatible with negligible optical depths in stars across a wide 
range of activity levels. Conditi ons in l a rge flar es may, how- 
ever, be different. Following iNess et alJ (|2003a), we measure 
line ratios between lines with low oscillator strengths, namely 
Fexvn /115.26 and Fexvn /116.78, and a line with a high os- 
cillator strength, namely Fexvn A 15.01. The flux in the lat- 
ter is thus more likely to be reduced by resonant scattering. 
Experiment ally, for zero opt ical depth /(15.26)//(15.01) = 
0.30 - 0.36 ( INess et alJl2003al and references therein). For the 
/(16.78)//(15.01) ratio, previous theoretical values range be- 
tween 0.40 - 0.50, while newer calculations bv lDoron & B ehar 
(2002) give a range between 0.60 - 0.75 depending on temper- 



ature (see also the detailed description in lNess et all2003al and 
references therein). 

The S/N ratio of our data was sufficient to derive some 
rough values for the above flux ratios. We fitted the lines with 
delta functions at their rest wavelengths convolved with the 
instrument response function. The continuum level was de- 
scribed, in the vicinity of the respective lines, by a power-law 
function. Since the line wings of the Fexvn /115.01 and the 
Fe xvn /115.26 overlap and an intervening Fe xix line is also de- 
tected at /115.20, these lines were simultaneously fitted. After a 
first iteration, the wavelengths of the brighter two lines (those 
at 15.01 A and 16.78A) were optimized by using them as fitting 
parameters as well. The deviations were minor. 

From our data, we find /( 1 5 .26)//( 1 5 .0 1 ) = 0. 1 8 - 0.3 1 for 
the intervals C, D, and E where the S/N is sufficient for mean- 
ingful ratios, with an error of about ±(0.06 - 0.09) for each 
ratio. For /(16.78)//(15.01), we find, for the same three inter- 
vals, values of 0.43-0.62, with errors of ±(0.06-0.12). We thus 
clearly do not detect any opacity effect at any of the flare phases 
investigated, at least within the large error ranges allowed by 
the present data. This suggests that even under extreme con- 
ditions of a flare plasma, optical depth effects are negligible, 
at least in the temperature range investigated with Fe xvn lines 
(T a 3 - 8 MK). 

5. Interpretation and discussion 

5.1. Small flares and coronal heating 

The first portion of our observation shows continuous variabil- 
ity to an extent that no part of the observation can be defined as 
being "quiescent" for intervals exceeding ss 1 hr (see Fig.0, 
and often variability occurs on much shorter time scales. This is 
a new aspect made possible by XMM-Newton 's h igh sensitivity 
and the long uninterrupted monitoring (see also lAudard et alJ 
2003b for a similar observation of the late-type M dwarf bi- 
nary UV Cet). We suggest that the numerous small peaks vis- 
ible in the X-ray light curve on time scales of 1-10 minutes 
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Fig. 8. Emission measure distributions for the pre-flare interval 
Q and for flare intervals A-E, using a polynomial reconstruc- 
tion algorithm that fits the DEM with polynomials of order 6. 
The lowest panel shows a superposition of all flare DEMs on 
the same absolute scale. 

are the most prominent examples of a large sequence of in- 
termediately luminous flares, approximately corresponding to 
M class flares in the solar classification. Support for this view 
is manifold: i) The light curves of the small enhancements re- 
semble typical flares on active stars or on the Sun (Paper I for 
examples), ii) We find many of these enhancements to be pre- 
ceded by optical bursts similar to the standard behavior in solar 
flares, and also similar to the correlated light curves of the large 
Proxima Centauri flare, iii) Hardness increases during periods 
of enhanced low-level activity. 

We suggested, in Paper I, that much of the X-ray emission 
observed during this observation may be due to the superposi- 
tion of a large number of intermediate and weak flares, most 
of which cannot be individually resolved, therefore forming a 



Fig. 9. Similar to Fig. [8] but using a regularization inversion of 
the spectrum to obtain the EM distribution. 

pseudo-continuous emission level with sporadic enhancements 
due to the strongest but least numerous members of the distri- 
bution. The principal model we thus propose is fundamentally 
different from a static corona and is rooted in evidence that 
magnetically active stars show mean coronal properties that are 
at variance with a corona like the Sun's. In a comparative study 
of solar analogs at different activity levels. lGiidel et al. (1997b) 
studied stellar X-ray EM distributions as a function of total X- 
ray luminosity (or age) and compared them with time-averaged 
EM distributions of solar flares. They found that active stellar 
EM distributions resemble the superposition of a low-level so- 
lar EM distribution and a time-averaged EM of strong (solar) 
flares as interpreted from GOES satellite data. Toward lower- 
activity stars, the hot component becomes weaker, probably 
due to a smaller rate of flares heating the corona to high temper- 
atures. The complete active-stellar EM distribut ions may thu s 
be the result of stochastic, superimposed flaring ( Gii dell 19971) . 
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In more active stars, the rate of strong flares is larger, but since 
larger flares heat more EM t o higher temperatures than do small 
flares (Feldm an et al the time-averaged EM distribution 

develops a prominent high-temperature component additional 
to the low-T emission measure. The latter is, in this picture, 
due to the combined effect of microflares that heat the plasma to 
only a few MK. The physical cause for an increase of the flare 
rate, and hence for an increase of the rate of large flares heating 
a large amount of plasma to high temperatures, was suggested 
to be th e increasing filling factor of active regions ( Giide l et alJ 
Il997bl) . As the concentration of active regions on an active star 
becomes larger, there are more interactions between tangled 
magnetic fields, leading to more numerous flares. The rate of 
very large flares thus also increases. 

This interpretation is also in line with recent studies of 
longer ti me series obtained with the EUVE an d BeppoSAX 
satellit es JAudard et all2000]lKashvap et al.l20ollGudel et"al] 
2003i see also lButler et alJll986l) that suggest that all of the 
observed extreme ultraviolet and X-ray emission of active late- 
type main-sequence stars could be made up by a population of 
stochastic flares distributed in total energy as a power law, with 
the smallest flares dominating the energy budget. All observed 
X-ray emission characteristics would thus be a consequence of 
a multitude of superimposed flare rises and decays. For exam- 
ple, the EM distribution would be dominated by the statistical 
properties (temperatur e, density, volume) a ttained during the 
long flare decays JGiidel et alJl997blEo03al) . 

Proxima Centauri is, during its low-level episodes, a mod- 
erately active star only, with Lx/Lboi ~ (0.6 - 2.4) x 1CT 4 , and 
with an EM distribution dominated by ^3 MK plasma and a 
tail reaching to rT O MK (Figs. EH 191. W e expect large flares 
to be much hotter. Feldman et alJ ( 119951) found a correlation 
between peak EM and peak temperature of solar and stellar 
flares. The faintest detectable flares in our observations were 
estimated to be (2 - 5) x 10 26 erg s -1 (Sect. 14. 1>. which im- 
plies EM * (7 - 16) x 10 48 cm 4 if a specific plasma emis- 
sivity per unit EM of 3 x 1CT 23 erg cm 3 s _1 is taken as a basis 
(from the SPEX software for temperatures of a few MK, us- 
ing solar abundances according to Table |3 for the flare peak; 
see lKaastra et alJ ll996a). Fig. 2 in Feldma n et alJ (11995) then 
suggests flare peak temperatures of 15 - 20 MK. Since these 
are only the largest of the contributing flares during low-level 
episodes and the characteristic temperatures during the long de- 
cay are typically a factor of two lower, our observations are 
well compatible with the view that low-level flares dominate 
the EM distribution. 

We can illuminate this problem from another angle. 
Extremely X-ray active M dwarfs with luminosities much 
higher than Proxima Centauri show higher characteristic coro- 
nal temperatures during low-level episodes. For example, the 
DEM of the active dMle+dMle binar y YY Gem is dominated 
by plasma at T = 6 - 12 MK JGiidel et alJl2001bh . The light 
curve is, somewhat similar to Proxima Centauri, interspersed 
with intermediate and weak flares, while a quasi-steady emis- 
sion level dominates the overall spectrum. The difference is 
that YY Gem's low-level luminosity is two orders of magni- 
tude higher. If the flare heating hypothesis as described above 
holds, then the low-level emission of YY Gem is dominated by 



0.002 



0.001 



Proxima Centauri (Flare) 



7 0.01 

E 





0.005 
0.004 
0.003 
0.002 

0.001 
0.000 



Proxima Centauri (Low) 



YY Gem 






10 15 20 25 

wavelength (A) 



30 



35 



Fig. 10. Comparison of the fluxed RGS spectrum of YY Gem 
(bottom panel) with the spectrum obtained during the large 
flare (combined intervals A, B, C, and E; middle panel) on 
Proxima Centauri, and the low-level spectrum (Q; top panel). 
The spectral bin width is 43.74A. 

proportionately larger flares, namely flares xTOO times stronger 
than those barely seen during Proxima Centauri's low-level 
emission, thus producing higher characteristic temperatures ac- 
cording to the Feldman et al. relation. Incidentally, the large 
flare on Proxima Centauri described here corresponds to such 
an event. We would then expect that the high-luminosity, long 
decay of the Proxima Centauri flare determines a spectrum that 
resembles the low-level emission of YY Gem. This is indeed 
the case. Fig.^|compares the fluxed RGS spectra of the com- 
bined flare portions A, B, C, and E (omitting the secondary 
flare peak) with the integrated sp ectrum of YY Gem collected 
from a * 1 d observation (Giidel et al. 2001b). The spectral 
line ratios and the continuum level look surprisingly similar 
but differ significantly from the much cooler spectrum observed 
during Proxima Centauri's low level emission. The flare DEM 
tend s to be so mewhat hotter than the integrated YY Gem DEM 
dGiidel et al.ll2001bl) . and the Neix and Nex lines at 13.45 A 
and at 12.1 A, respectively, tend to be somewhat lower in the 
Proxima Centauri flare compared to the Fexvn lines, which 
is probably due to differences in the Ne/Fe abundance ratio 
(Fexvn and Neix have very similar line formation tempera- 
tures ). Note, however, t hat the DEM reconstruction presented 
bv lGiidel et all (12001 bl) was based on RGS2 only; this spec- 
trum is rather insensitive to temperatures significantly above 
10 MK. The MOS spectra (Table 2 in lGudel et al1l2001bl) in- 
dicate an additional large amount of EM around 20 MK (see 
also Stel zer et al.l 12002). The presence of weaker flares and 
later flare decays will also increase somewhat cooler compo- 
nents in the DEM. This comparison again suggests that more 
active stars are dominated by flare-like plasma correspond- 
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ing to intermediate-to-large flares, reaching to high tempera- 
tures, while lower-activity coronal emission may as well be 
composed of flare contributions which, given the lower flare 
energies, attain lower average temperatures. This hypothesis 
would explain i) why more active stars show hotter coronae, 
ii) why more active stars show much higher luminosities (due 
to the increased EM built up during the flare process). It fur- 
ther supports the model that active stellar coronae are heated by 
stochastic flares. And lastly, it may explain why active (main- 
sequence) stars show high characteristic densities. If low-level 
emission is made up of flares as the one discussed in this pa- 
per, then densities derived spectroscopically from the time- 
integrated spectrum across the entire flare should be similar to 
densities measured in active stars during quiescence. For the 
Proxima Centauri flare, we obtained an average characteris- 
tic density from the O vn triplet of log« e = 10.50+" while, 
for example, the low-level YY Gem spectrum reveals a simi- 
lar value, namely lo gw,, _ 10.35+2 jj, as do many other a ctive 
main-sequence stars (Gi idel et all2001allNess et alJ2002h . 

5.2. Flare models 

We will now discuss the large flare in the context of standard 
flare models. Various models have been discussed in the liter- 
ature, although all of them represent simplified paradigms that 
may not be realized in the Sun or stars in their pure forms. 
Rather, we expect that some gross characteristics of the models 
may approximately describe the observations, thus providing 
some insight into the overall development of magnetic energy 
release. With this in mind, we will, in the following, describe 
two complementary models and investigate their applicability 
to large flares as the one presented her e. The first approach, a 
2-ribbon flare model based on work by Kopp &Polettol ( fl98l 
but extended to include approximate radiative and conduc- 
tive c ooling losses, was presented in this form bv lGiidel et alJ 
The second approach, using full hydrodynamic simu- 
lations of a closed loop or a pair of l oops, has been previously 
described by, e.g. JPeres et alJ I I1982I) . We will discuss one typ- 
ical numerical solution, but will defer an in-depth presenta- 
tion of the hydro dynamic simulations to a companion paper 
jReale et all200 4). We begin this discussion by briefly review- 
ing some information on the initial energy release phase that 
can be extracted from the optical burst. 

5.2.1 . The impulsive phase of the flares 

The initial energy release of a solar flare is in general most ob- 
viously seen at radio wavelengths, in the optical/UV regime, 
and in non-thermal hard X-rays. The relevant signals in these 
three wavelength regimes are thought to be related to high- 
energy electrons accelerated in the initial phase of a flare in the 
coronal reconnection region. While radio signals are due to gy- 
rosynchrotron emission of electrons spiralling in the magnetic 
fields, the hard X-ray signal is due to bremsstrahlung when 
a non-thermal population of electrons impacts in the chro- 
mospheric/transition region layers at the magnetic field foot- 
points. At the same location, optical emission is thought to 



arise promptly due to heating and ionization of the chromo- 
spheric plasma. As this outline suggests, the three signatures 
of the impulsive flare phase should occur early in the flare, as 
the electrons are thought to be the cause for the subsequent 
chromosp heric evaporation that then dominates the soft X-ray 
emission jAntonucci et al.l["l984l |Pennislli988l) . For the opti- 
cal emission observed in Proxima Centauri, this is clearly the 
case as was shown in Paper I (see also Fig.[2j. Furthermore, so- 
lar observations clearly show - as expected - a very close time 
correlati on between the hard X-rays and the microwave emis- 
sion ( iKosugi et al . 1988), and between optical "white light" and 
hard X-ray emission faudson et alJl992HDennisll98 8 and ref- 
erences therein), often within seconds. 

While solar research has made extensive use of non-thermal 
hard X -rays for the study of the energetics of the impulsive flare 
phase (Dennis 1988), the respective signals are presently out of 
reach for stellar astronomy since required detector sensitivities 
are not available at the energies of interest (>20 keV). However, 
observation of "white light" flares (or equivalently, as in this 
work) U band flares is relatively easy for M dwarfs given the 
large contrast of flares against photospheric background emis- 
sion. The U band has traditionally played a major role since 
the flare enhancement is particularly large there. The physics 
of the continuum enhancement is not very well understood. 
iHudson et alJ i 19921) (and references therein) argue for exci- 
tation by energy transfer from electron beams, either by direct 
collisions, or by UV irradiation excited by the el ectron impact. 
For a further, short discussion on the issue, see IHawlev et alJ 
d!995h . 

Independent of the poorly understood physical cause of the 
U band emission, our observation, together with previous so- 
lar knowledge puts us into a position to at least grossly assess 
the time evolution of the primary energy release in electrons 
if the standard model outlined above holds. The U-band light 
curve (above its photospheric level) would thus be a sensitive 
tracer of the production of non-thermal particles in the corona, 
and therefore an efficient proxy of the inaccessible non-thermal 
hard X-ray emission so often seen in solar flares. 

A second application is possible, namely to roughly de- 
termine the chromospher ic footpoint area of the flaring mag- 
netic fields, as outlined in lHawlev et al.Nl995l) . To do this, one 
would have to know the temperature of the U-band emitting 
source, requiring at least one more observation i n another opti- 
cal ban d, which is not available. The analysis by IHawlev et alJ 
d 1995b for flares on the M dwarf AD Leo suggests tempera- 
tures of as 9000 - 9500 K under the assumption of blackbody 
emission. Assume that the filling factor of the flare footpoints 
in the projection onto the sky is X. Since the area emitting at the 
non-flaring photospheric temperature is reduced by the flaring 
area, the U-band luminosity ratio between measurements taken 
during the flare and outside the flare is 

D L f + (l-X)L q XB x] (T f ) + (\-X)B v (T q ) 

R=— - = — . (1) 

Lq B\]{T q ) 

Here, Lf is the U-band flare luminosity, and L q is the 

U band quiescent luminosity of the whole star; B\j is the 
Planck function evaluated in the U band 3500A) for 

the flare temperature (Tf = 9000 K, see above) and for 
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the photospheric effective temperature outside flares (T q = 
2700 K. iFrogel et al.lll972l) . respectively. Therefore, recalling 
that B V (T) oc (exp(hc/kAT) - 1)~ 1 « exp(-41 156/J), we find 

R = 1 + x(e 4U56(1/T «- l/T ^ - l) * 1 + 4.3 x 10 4 X. (2) 

Our flare increased the U band flux by a factor of R a 20, so 
that the fractional projected area of the footpoints is « 4 x 10~ 4 
of the cross section of the star. If this described one circular 
footpoint, then its diameter would be 2% of the stellar diam- 
eter. Comparing with the results we will derive from X-rays 
below, the area covered by U-band emitting footpoints must 
be considerably smaller than the X-ray emitting area, which is 
compatible with magnetic loops that span over a wide range 
across the stellar surface. 

5.2.2. A two-ribbon flare model 

The two-ribbon (2-R) flare model devised initially by 
iKopp & Polettol Jl984h is intrinsically a parameterized 
magnetic-energy release model in which the time development 
of the flare light curve is fully determined by the amount of 
energy available in non-potential magnetic fields, and the rate 
of energy release as a function of time and geometry. The 
complete time evolution of the flare is based on continuous 
heating that feeds an evolving reservoir of hot plasma. Plasma 
cooling is not included in the original model; it is assumed 
that a portion of the total energy is radiated into the observed 
X-ray band, while the remaining energy will be lost by other 
mechanisms. 2-R flares are well established for the Sun; they 
often lead to large, long duration flares that may be accompa- 
nied by mass ejections. This configuration is clearly suggestive 
for our observation, paralleling conclusio ns for similar flare s 
observed previously on Proxima Centauri (Hais ch et alJl9 83). 
We will describe here an extension of the model that includes 
approximations of radiative and conductive losses, and that 
will also take into account some geometric development of the 
radiating magnetic structures. This variant of the 2-R model 
was first applied to a large, lon g-duration flare ob served on the 
RS CVn-type binary UX Ari jGudel et all 1 1999b . A detailed 
description is given in the latter reference; we give a brief 
summary below. 

A 2-R flare conceptually starts when a disruptive event, 
e.g., a filament eruption, opens a loop arcade. The open field 
lines are then driven toward a radial neutral sheet (above the 
magnetic neutral line) where they reconnect at progressively 
higher altitudes, forming a system of overlying loops below 
them. The excess energy provided by annihilation of the mag- 
netic fields is released in an unspecified form. The amount 
of continuous heating provided throughout the reconnection 
episode principally describes both the flare increase and its de- 
cay. 

The magnetic fields are, for convenience, described along 
meridional planes on the star by Legendre polynomials P n of 
order n, up to the height of the neutral point; above this level, 
the B field is directed radially. The magnetic fields are thus as- 
sumed to be poloidal, initially directed radially outward and 
then reconnecting, starting with the antiparallel field lines clos- 
est to the neutral line. As time proceeds, field lines further away 



from the neutral line move inward at coronal levels and recon- 
nect at progressively larger heights above the neutral line. The 
reconnection point thus moves upward as the flare proceeds, 
leaving closed magnetic loop systems underneath. One loop 
arcade thus corresponds to one (N-S aligned) lobe between 
two zeros of P„ in latitude, axisymmetrically continued over 
some longitude in E-W direction. The propagation of the neu- 
tral point in height, y(t), with a time constant frj (in units of R*, 
measured from the star's center) is prescribed by 



y(t) = 1 + 



R. 



H(t) = [y(t)-l]R t 



(3) 



(4) 



and the total energy release of the reconnecting arcade per ra- 
dian in longitude is equal to the magnetic energy lost by recon- 
nection, 

= 2n(n + l)(2n + l) 2 /?fe 2 / 12 ( n ) r y ) y ' (5) 

dy 87T [n + (n+ l)y 2n+l ] 3 



dE 
~dt~ 



dEdy 
dy dt 



(6) 



(Poletto et al. 1988). In Eq. H m is the maximum height of 
the neutral point for t — > oo; typically, H m is assumed to be 
equal to the latitudinal extent of the loops, i.e., 



n+l/2 



-R, 



(7) 



for n > 2 and H m = (7r/2)R t for n = 2. Here, B is 
the surface magnetic field strength at the axis of symme- 
try, and R„ is the stellar radius. Finally, Inin) corresponds to 
J [P n (cosff)] 2 d(cos8) evaluated between the latitudinal borders 
of the lobe (zeros of dP„/dd), and 9 is the co-latitude. We have 
evaluated this expression for lobes that are centered at the star's 
equator for odd n, and lobes that start at the equator for even n. 
Note that R*, B 2 , and Iniri) merely define the normalization of 
the energy release light curve, not its form. The latter is deter- 
mined by the time constant to and the polynomial order n. 

The largest realistic two-ribbon flare model is based on the 
Legendre polynomial of degree n — 2; the loop arcade then 
stretches out between the equator and the stellar poles. For sim- 
plicity, we assume that the arcade extends along a great circle 
on the stellar surface, e.g., the equator (instead of adopting ax- 
isymmetric arcades at any other stellar latitude), although the 
computations of the latitudi nal and radial geome try use the ex- 
act equations as given by IKopp & Polettol dl984). Since solar 
two-ribbon flares occur in loop arcad es whose leng th is typi- 
cally about 1.5 times their width (Pole tto et alJl988l) . we adopt 
this ratio for any given n: L — l.5H m = 3R t n/(2n + 1). 

The Kopp & Poletto model is applicable after the initial 
flare t rigger mechanism has terminated, although IPneumar] 
dl982l) suggested that reconnection may start in the earliest 
phase of loop structure development. The time of strong pro- 
duction of non-thermal electrons, signified by the optical burst 
during the X-ray increase, gives rise to explosive chromo- 
spheric evaporation and steeply increasing temperatures. This 
phase of the flare may be less well described by the model. 
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IPoletto et all dl988h introduce an ad hoc factor of 10% of 
the magnetically released radiative energy that escapes into the 
X-ray regime. Here, we include time-dependent conductive and 
radiative losses in the X-ray corona self-consistently and hence 
account for the emitting volume. The volume and the observed 
EM provide an overall electron density for which we obtained 
estimates from the He-like line triplets in the RGS spectra. 
Requiring that the densities agree with the measured densities 
will thus constrain the model parameters, in particular n. The 
radiative loss time is 



r r = 



3kT 
n e ip 



(8) 



where if/ is the emissivity per unit EM of the hot, optically thin 
plasma, which we take to be if/ = 2.5 x 10~ 23 ergs s^'cm 3 for 
T = 15-20 MK (from the collisi onal ionization equilibrium 
model in the SPEX software, see iKaastra et alJll996ah . The 
downward conductive flux in a constant cross section loop is 



ds 



dT 



(9) 



with s being the length parameter along the loop, and k a 
10~ 6 ergs cirr's 'K -7 / 2 dSpitzerl 1 1 962h . F c can be approxi- 
mated as 



n7/2 



F c ^ TjK- 



(10) 



where th e peak temperatu r e T m » T, tj — 4/7 has been intro- 
duced hv lKonn&Polettol {1993 ) to account for the geometry 
in this approximation, and I is the time-dependent semi-length 
of the loop, 



/(f) * -Hit). 



(ID 



Cooling is principally controlled by radiation and conduction 
with time constants of r r and t c , respectively, leading to an 
effective cooling time r e given by 



1 _ 1 1 

T e T c T r 



(12) 



Combining the relevant equations (see Gii del et alJl999l for de- 
tails) one obtains for the model X-ray luminosity 



dE L 

L x (t) = f(t) 

and the underlying EM must be 



EM(f) = f(t) 



dE/dt L 
W) R~* 



(13) 



(14) 



Here, fit) = r e (f)/r r (f) gives the fraction of the total energy 
being emitted as X-rays, and L/R* denotes the length of the ar- 
cade in radians, here taken to be < 2 (L = 1 .5H m ). We adopt T 
from the observations, T » 20 MK during the period of consid- 
eration. Note that the observable X-ray emission is not propor- 
tional to dE/dt in Eq. (1131 as the fraction / is time dependent. 
The run of the electron density n e then is 



«e(0 = 



dE/dy 
3kTR t nHit) ' 



(15) 



The electron density is solely determined by the geometry of 
the arcade, as prescribed by Eq. Q-Q. 

We integrated all relevant equations in short time steps for 
a given n, adjusting to and B until a good fit to the observa- 
tions was obtained. Following Giidel et al. ( 1999), we assume 
that only a fraction q < 1 of the liberated magnetic energy is 
used to heat the plasma (energy that is lost via radiation and 
conduction, while the remaining fraction is transformed into 
mechanic al energy, into fast particles ejected from the corona, 
etc). Kopp & Poletto ( 1984) use q — 0.003 in their application 
to a solar flare. 

As the overlying magnetic field lines reconnect, the loops 
already reconnected are pushed downward until they asymptot- 
ically reach a final position, corresponding to a potential field 
configuration. We have assumed that the X-ray loops have re- 
laxed to this final geometry while radiating most of their X-ray 
energy (corresponding to t — > oo). The derivation of the asymp- 
totic loop conf iguration is p erformed by piecewise integration 
as described in Giid el et alJ (11999). 

There are only few independent parameters for this model. 
The total energy released is controlled by B 2 , and by the effi- 
ciency q. To obtain the same thermal energy, q obviously scales 
as q oc B 2 . The fraction of radiated X-ray energy out of the to- 
tal thermal energy is given by fit). These factors together thus 
determine the X-ray flare amplitude. The time scale is given by 
the parameter to in Eq. l|3}. Additionally, Eq. dl5> provides an 
approximation of the loop density that can be compared with 
the observations. There is thus a family of solutions (B, q) for 
each n. The best-fit solution for n — 2 is shown in Fig. \^\ 
(top) and required a maximum surface magnetic field strength 
of B - 2050 G for q = 0.01. A larger efficiency of q = 0.05 
thus corresponds to B = 920 G. These two examples charac- 
terize the probable range for B although stron ger fields have 
also b een measured in late-type dwarfs Ijohns-Krull & Valentil 
2000). The time constant f ~ 4000 s, and the final height of 
the neutral point (for t -» oo) is H m « 1.73 x 10 1() cm or 1.71?, 
for n - 2. 

For n = 3, we find B = 2720 G for q = 0.01, and t = 
3300 s. Similarly, for n = 4, the best-fit parameters are, for 
q = 0.02: B = 2200 G and f = 3100 s. Since we can also 
estimate the model densities from Eq. dl5> . our observations 
of density-sensitive spectral lines, in particular O vn, provide a 
further constraint. For n = 2, Fig.[^ s h ows an electron density 
at flare peak of n e « 3 x 10 11 cm -3 , quite compatible with 
the O vn measurements. The latter refer to much cooler plasma 
than the bulk plasma in a flaring loop at its peak luminosity. 
However, as discussed in Sect. 14.31 densities inferred for hot 
solar flare plasma seem to be of similar magnitude as densities 
derived from the Ovn lines. For n — 3, densities at flare peak 
are « 5 X 10 11 cirT 3 , while for n = 4, n e « 7 x 10 11 cirT 3 . In 
the latter case, the magnetic field estimated from the model at 
the loop top and at the time of flare peak is below the pressure 
equilibrium value if the electron temperature is « 20 MK. We 
therefore tend to exclude solutions for n — 4 and higher within 
the context of this simple model. 

The model fits well until about t = 1000 s after flare start; 
thereafter, the model light curve clearly decays too fast, which 
is due to the emergence of the secondary flare. Indeed, at t » 
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pared with the evolution of the flare. If this were not the case, 
then a significant portion of the energy released at earlier times 
would contribute to the emission at time t, i.e., the light curve 
would be a convolution with radiation from previous times. 
From Fig. |2] we estimate that changes in the light curve oc- 
cur on a time scale of about 1000 s near flare peak. The fit is 
satisfactory if r e r 1000 s up to the flare peak. Our evaluation 
of Eqn. H2\ indicates that r e = 80 - 100 s around the flare peak 
for the acceptable solution presented in Fig. ^2 an d the radia- 
tive loss time r r is around 800-1000 s, which is compatible 
with our requirements. 

We caution that the assumptions for this flare model are 
crude. The purpose of our estimates, however, is to indicate 
whether the observed flare belongs to a class of compact flares 
or may be a global phenomenon, compatible with the standard 
two-ribbon flare model. Our results thus suggest that the flare is 
likely to be geometrically large, with size scales (height, foot- 
point separation) on the order of 10 10 cm or one stellar radius. 
The presence of such extended coronal structures is compati- 
ble with VLBI observations that show flaring cores that rapidly 
evolve into extended structures dMutel et alJll985|) . Observed 
radio so urces reveal size scales on the order of the stellar di- 
ameters (Ben z et alJll998l) . Also, the absence of radio eclipses 
or strong signatures of rotational modulation in long-duration 
flares on similar M dwarfs agrees well with the assumpt ion of 
extended flare geometries on such stars jAlef et all 19 97). 



Fig. 11. Top: Simulation of the large Proxima Centauri flare 
using a simple 2-R flare model described in the text. The best 
fit shown here is for n = 2, B = 2050, f = 4000, q = 0.01. 
The solid line shows the model X-ray luminosity, to be com- 
pared with the observations that are plotted by small crosses. 
The dotted curve indicates the electron density (see scale on 
the right), and the two solid lines at the bottom of the figure 
mark the height of the emitting loop (see scale on the right). 
The model clearly decays too fast beyond t = 1000 s after 
flare start, which is due to the emergence of a secondary flare. 
- Bottom: Similar to upper figure, but here the light curve syn- 
thesized from a time-dependent hydrodynamic simulation of 
plasma confined in a single footpoint-heated loop is overplot- 
ted (solid line). The heating lasts 600 s. The flare start is defined 
slightly differently to optimize the fit to the rise phase (the pre- 
cise start of the flare cannot be defined from the soft X-rays 
unam biguously). Further details are described in iReale et all 
(2004)/ 

1000 s, the hardness ceases to decay and dev elops a plateau (se e 
Fig. |3, presumably due to resumed heating jReale et alJ20 04). 
We emphasize that in the framework of this model, the flare 
density is not a consequence of cooling material condensing 
out of the corona and falling back to the surface. Rather, the 
loops radiating at later times are less dense than those radiating 
at earlier times, according to Eq. J15I . The 2-R model does not 
describe the evolution of plasma within a loop, but rather the 
plasma and geometry parameters of a series of loops in time. 

We have to check the value of the decay time r e . The model 
presented here assumes that the decay time r e is short com- 



5.2.3. Hydrodynamic flaring loop model 

X-ray flares can alternatively be modeled using full numeri- 
cal hydrodynamic simulations. We outli ne the meth od here, but 
for a technical des cription, we refer to Peres et all dl982t) and 
IReale etalJ (Eo04). The latter reference discusses applications 
to the present observation in deeper detail. The numerical ap- 
proach stresses on plasma confinement. The basic assumption 
becomes that the flare is mostly governed by the evolution of 
the plasma confined in a magnetic loop and that changes of the 
loop system geometry and size are small, slow and play a minor 
role. 

One can then focus on modeling the flaring plasma con- 
fined in a single unchanging loop, where it moves and and 
transports energy along the magnetic field lines. The magnetic 
field does not play any role directly, except for plasma confine- 
ment. 

The confined plasma behaves as a compressible fluid and 
its evolution can be described with a one-dimensional time- 
dependent hydrodynamic model. When solving the equations 
of conservation of mass, momentum and energy, in coronal 
plasma conditions, several physical terms are important, and, in 
particular, the pressure gradients, the stellar gravity, the radia- 
tive losses and the thermal conduction. Describing all of them 
requires the numerical solution of the hydrodynamic equations. 
A parameterized and time-dependent heating function is in- 
cluded to trigger the flare. A sufficiently intense heat pulse is 
able to make the plasma temperature increase to flare values 
and to drive a strong plasma evaporation from the chromo- 
sphere: the combination of these two effects switches the flare 
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on. The flare decay is instead mostly driven by the plasma cool- 
ing by radiation and thermal conduction, after the heating pulse 
has finished, unless a residual heating is able to sustain the de- 
cay. This makes an important conceptual difference from stan- 
dard 2-R models, in which the decay is assumed to be totally 
driven by the release of magnetic energy during the reconnec- 
tion of progressively larger arches. 

Several numerical time-dependent loop models have been 
applied to study solar coronal flares. One of them has been ap- 
plied to study part of a flare observed on Proxima Centauri with 
the Einstein satellite jReale etalJl988h . 

In an accompanying paper we model the flare for most of 
its duration, including the second peak, by describing the evo- 
lution of flaring plasma confined in single loop structures. The 
detailed numerical approach allows us to synthesize the plasma 
X-ray emission at the focal plane of the EPIC PN and to make 
detailed comparisons with the real data. We present the overall 
results of a simulation here for comparison with our analytic 
approach using the 2-R formulation. 

In the simulations, the flaring plasma confined in a single 
loop is able to explain well the rise phase, main peak and initial 
decay of the flare (Fig.[^ bottom), and the addition of an ar- 
cade of similar flaring loops allows to explain the second peak 
and the late decay. In spite of the different approach, the results 
are not in contrast to those obtained with the 2-R modeling 
described in Sect. 15.2.21 Hydrodynamic modeling and detailed 
comparison of results to data will however allow us to obtain 
specific constraints on the loop geometry and on the heating 
function throughout the flare, and insight in the physical flaring 
plasma conditions. In particular, the loop height appropriate for 
the simulations is very similar to the typical height found in the 
2-R model, namely tss 0.6 x 10 10 cm. This result is thus in good 
agreement with estimates based on the 2-R model. The princi- 
pal finding from the modeling is thus that the flaring source is 
intrinsically large, of order of one stellar radius. 

6. Conclusions 

Proxima Centauri is the closest star to the Sun, and it belongs 
to the important class of magnetically active stars as well. The 
latter are generally characterized by enhanced spot coverage in- 
dicating a large surface magnetic filling factor, chromospheric 
emission lines such as Ha, enhanced flaring activity, or strong 
coronal emissions. Proxima Centauri does not show signatures 
of extreme activity but given its surface area of about 2% the 
Sun's, its X-ray luminosity equaling an average solar level is 
notable. A more important aspect of the present and all previ- 
ous X-ray investigations of this star is in fact simply its very 
late spectral type, which has two important implications : First, 
since optical flare emission is predominantly "blue", a large 
contrast against the weak, red photospheric emission makes op- 
tical flare detection easy. Our OM data have been an important 
guide to recognize low-level flaring in the X-ray light curve. 
And second, there is quite limited surface area to host active 
regions of a size that would be typical for the Sun. In fact, 
Proxima Centauri is no larger than a large solar coronal active 
region JReale et all20 04). This circumstance makes X-ray ob- 
servations very sensitive to reconfigurations in the corona that 



likely affect a considerable portion of the surface. This obvi- 
ously helped set a good contrast for the large number of weak 
flares we have detected during the low-level X-ray emission. 

Our observation has captured a continuously variable X-ray 
corona, with flares spanning two orders of magnitude in peak 
flux and at least three orders of magnitude in total X-ray en- 
ergy. The weakest flares seen during the low-level episode are 
principally the weakest of their kind observed on any active 
star, while the large flare is the largest yet seen on Proxima 
Centauri, exceeding even some of the strongest solar X-ray 
flares. Despite this large range of flare parameters, we find 
similar timing characteristics among them: Many X-ray flares 
are preceded by optical bursts, suggesting the process of chro- 
mospheric evaporation by electron bombardment as the main 
driver of heating and mass transport into the corona. The large 
flare is only set apart by its very slow decay and secondary 
peaks that indicate continued heating. 

We interpret the low-level emission as being predominantly 
due to a superposition of stochastic flaring. While a truly quies- 
cent emission level due to steadily heated coronal loops cannot 
be excluded, the continuous decay of emission by a factor of 
ten during the first 12 hours of our observation requires rela- 
tively rapid reconfiguration of the corona unless the complete 
flux decay is due to the late decay of a large flare that reached 
its peak before our observation started. (Rotational modulation 
of active regions is not likely to be important on the observed 
time scales; the rotation period of Proxima Centauri i s sus- 
pected to be as long as 83.5 d, see B enedict et alJll99 8.) The 
additional modulation on time scales of hours and associated 
hardness changes argue for physical changes in the corona on 
such time scales. If most of the low-level emission originates 
from flaring plasma, then the flare peaks individually detected 
constitute the strongest examples of a population of low-level 
flares. As such, they determine the high-temperature cutoff of 
the emission measure distribution since flare peak temperatures 
are d ependent on the peak emission measure JFeldman et alJ 
1 19951) . We have argued that the expected temperatures for these 
flares are roughly commensurate with the observed EM distri- 
bution: the DEM shows a high-T extension to » 15 MK, which 
is the expected peak temperature for the largest flare during the 
low-level emission. 

If this picture has some merit, then it should apply to other 
stars at higher activity levels. For a characteristic distribution of 
flare energies as seen on the Sun, higher flare rates also induce 
larger rates of larger flares. Therefore, a larger flare rate (above 
a certain threshold energy) brings more material into the corona 
and at the same time heats some of it to much higher temper- 
atures according to the Feldman et al. relation. We thus expect 
characteristic temperatures in more active stars to be higher, 
which is generally the case. We have scaled this model to the 
extremely active dMe binary YY Gem, which shows an X-ray 
luminosity »T00-200 times higher than Proxima Centauri. The 
largest flares frequently dominating the low-level emission of 
YY Gem are consequently similar in energy to the large flare 
in Proxima Centauri that we captured during the second half of 
our observation. The combined spectrum of the latter flare in- 
deed resembles the average spectrum of Y Y Gem, and so do the 
EM distributions although a full modeling should take into ac- 
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count the additional lower-energy flares occurring on YY Gem 
as well. 

We therefore suggest that the low-level emission of mag- 
netically active stars is made up, to a large fraction, by a su- 
perposition of flares over a large range of total energies. The 
flare rate (above a given low-energy threshold) determines the 
high-temperature end of the EM distribution of the intervals 
perceived as "low-level". A natural consequence of this model 
is that, as the flare rate increases toward more active stars, the 
emission measure increases because more chromospheric evap- 
oration is taking place, while the large amount of EM produced 
by larger flares is also heated to higher temperatures, thus in- 
ducing a hotter EM distribution in these stars. We thus suggest 
that the Feld man et relation naturally explains why 

stars at a higher activity level show hotter coronae. 

The large flare observed during the second part of our ob- 
servations has offered a number of new views on extremely 
large stellar flares. Like most of the other, low-level flares, the 
earliest phase of the rapidly heating X-ray episode is accompa- 
nied by a very large optical burst. From solar analogy, we in- 
terpret this "white light" flare as being the signature of acceler- 
ated electrons bombarding the chromosphere and evaporating 
the material that is subsequently becoming visible in X-rays. 
In this sense, we find that U band flares are ideal and sensi- 
tive proxies for the non-thermal hard X-rays (> 20 keV) often 
investigated in solar flare research. 

While the emission mechanism of the optical flare is some- 
what uncertain, we infer relatively small areas for the chromo- 
spheric footpoints. This is in contrast to the X-ray geometry we 
find from two different approaches (an analytic 2-R model and 
a fully hydrodynamic approach). Both X-ray flare models indi- 
cate relatively large flares, with typical size scales of 10 10 cm 
or one stellar radius. While Proxima Centauri is much smaller 
than the Sun, the size of such active regions is comparable to 
those that produce very large flares on the Sun. 

For the first time, X-ray spectroscopic measurements have 
fully resolved the electron density development during a stellar 
flare, indicating peak densities in the 1-5 MK plasma of ap- 
proximately 4 x 10 11 cirT 3 . Such values are similar to densities 
measured during solar flares. We have also measured the abun- 
dance development during several stages of the flare. Apart 
from a trend toward higher abundances of high-FIP elements, 
the abundance ratios are relatively close to solar photospheric 
ratios. 
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